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O ■ Abstract 

o\ : 

We investigate hard exclusive reactions on the nucleon with soft pion emission. A 
parametrization of corresponding hadronic matrix elements in terms of parton distri- 
butions for final pion-nucleon state is provided. These distributions are calculated in 
terms of nucleon and pion GPDs and the pion distribution amplitude via soft-pion 
theorems. Some observables for the process of hard charged pion production on the 
proton with soft pion emission are computed. 
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1 Motivation and outline 

The field of hard exclusive reactions has been studied intensively during the past decade. 
The investigations have proceeded in the theoretical as well as in the experimental sector; 
reviews are given e.g. in references [H El EH IU E| • Two prominent representatives of the 
hard exclusive processes are deeply virtual Compton scattering (DVCS) and hard meson 
production (HMP). The nucleon properties which enter these reactions are formulated in 
terms of generalized parton distributions (GPDs). On the one hand, these functions can 
be viewed as generalizations of the usual forward parton distributions, on the other hand 
they are directly related to nucleon form factors through their moments. So it has been 
argued that hard exclusive reactions can provide useful new insights into the partonic 
nucleon structure which are not accessible through the usual electroweak probes. In this 
context, for example, the form factors of the energy-momentum tensor have been discussed, 
see e.g. references Additionally, in the case of meson production, one can obtain 

information about the involved distributions amplitudes. 

In this article, we investigate the situation when in a hard exclusive reaction instead of 
the final nucleon a nucleon-pion state with low invariant mass appears. Since it is produced 
close to the threshold, the pion is denoted as soft. On the experimental side, a separation of 
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a particular hard reaction with and without soft pion emission cannot always be guaranteed. 
In this sense, the process with soft pion can be viewed as a contamination of the fully 
exclusive DVCS or HMP, and theoretical estimates about this disturbing background are 
desirable. 

Apart from such practical considerations, these new reactions are worth being studied 
in their own right. They provide an opportunity to investigate soft pion emission from the 
nucleon induced by nonlocal lightcone operators as opposed to the local vector or axial 
operators to which we are restricted in usual electroweak pion production. Therefore, 
analogously to how a soft process such as pion-electroproduction can provide information 
about nucleon form factors, hard processes with soft pion emission might contribute to a 
better understanding of quantities such as generalized parton distributions. 

Guichon et. al. have addressed this question for the process of DVCS with soft pion 
production |Hj. For the calculation of the pertinent hadronic matrix element, they presented 
a soft-pion theorem based on current algebra and chiral symmetry. Moreover, they modeled 
the effect of the A(1232) resonance which is located not far from the pion production 
threshold. By this method, they could give predictions for certain cross sections and 
asymmetries. However, in all their considerations, the region of small momentum transfer 
was explicitely excluded. 

In the following, we present our approach of calculating soft pion emission in hard 
exclusive reactions. Based on polology arguments, "PCAC", and certain properties of the 
chiral symmetry transformation, we derive corresponding soft-pion theorems. These differ 
from the results in |H] through additional pion pole terms. The implementation of these 
new contributions allows in particular to extend the region of applicability down to small 
momentum transfer. Using this improved result, we calculate the effect of soft pions in 
hard 7r + production off the proton. 

The outline is as follows. Basic kinematical considerations are given in section 2. In 
section 3, we provide a parametrization of the matrix elements for pion emission induced 
by twist-2 lightcone operators. We denote the invariant functions that come up in this pro- 
cedure as pion-nucleon (nN) parton distributions. These functions are the generalizations 
of GPDs for the case of pion emission. We discuss some of their properties, in particular 
their behavior at the pion threshold and the meaning of some moments. In section 4, we 
give a detailed derivation and the results of soft-pion theorems for several twist-2 operators. 
We check that in certain limiting cases our expressions are consistent with previous calcu- 
lations. Finally, in section 5, we apply these results to the process of hard pion production 
off the proton with soft pion emission near threshold. The amplitude of the process for 
arbitrary isospins is given. Further, we provide the following numerical estimates for hard 
7r + production: the transverse spin asymmetry of the process 7* + p — > ii + + N + 7r S oft is 
calculated, and the contamination of the longitudinal cross section and the transverse spin 
asymmetry of the pure process 7* + p — > vr + + n through soft-pion admixture is determined. 
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N(p,S)' t N(p'.S') 



Figure 1: Particle momenta and invariants for the exclusive reaction with soft pion emis- 
sion. 



2 Kinematics 

We consider the collision of a virtual photon 7* with momentum q and a nucleon N with 
momentum p and spin 5*. In the final state, we have a nucleon with momentum p' and 
spin S', a pion tt with momentum k and isospin a, and either a real photon 7 (DVCS) or 
a specified meson M with momentum q'\ 



1 *( q )+N(p,S)^N(p>,S') + 7r a (k) + 



!(</) 
M(q') ' 



It is useful to define an average momentum p and a momentum transfer A in the following 
way: 

_ p + p' + k , , 

p = , A = q-q=p+k-p. (1) 

Further, we introduce the Lorentz invariants 

s = {p + q)\ t = {p-p')\ u = (p-k) 2 , W 2 ^{p' + k)\ (2) 

see also the illustration in figure [T] (we point to the fact that here A 2 is not identical to 
the nucleon momentum transfer t unlike in the case without pion). The kinematical region 
of hard scattering is characterized by a large photon virtuality Q 2 , a large energy v of the 
virtual photon in the target rest frame, and a fixed Bjorken variable xb- 

Q 2 = -q 2 , Mv=p-q > -A 2 , W 2 , M 2 , m 2 , q' 2 , (3) 
Here, M and m, denote the nucleon and pion mass, respectively. 
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We shall refer to the pion ir(k) as soft if it appears sufficiently close to its production 
threshold in the following sense. Directly at the threshold, the variables W 2 and u are 
fixed as 

(5) 



Wl = (M + m^) 2 , u th — (M — m n ) 2 + -^t. 



We allow for deviations from these values that are of the size 

W 2 - W 2 h < Mm n , (u - u th ) < Mm n . 

Moreover, we require for the nucleon momentum transfer that 

-t < M 2 . 



(6) 



(7) 



Relations @ to settle the soft pion kinematics (note that they imply —A 2 < M 2 ). 
Roughly speaking, we can summarize these conditions simply as 



k = 0(m„), p,p' = 0(M). 



(8) 



Let us now turn to a Sudakov decomposition of the particle momenta. Two lightcone 
vectors n and n are defined in a frame where p = (po, 0, 0, p z ) and q = (go, 0, 0, q z ) via 



n= (Po +&)(!, 0,0,1) 



n 



2(Po + Pz 



;i, 0,0,-1). 



(9) 



Consequently, we obtain the average momentum and the virtual photon momentum in the 
form 

M 2 Q 2 
p = n + —n, q = -2£n + — n, (10) 



with 



-, 2 M 2 + W 2 A 2 

= p 2 T' 



-2e 



P-Q 

M 2 



f P-Q 

\M 2 



\ Q 2 



M 2 



The decompositions of some other vectors to leading order in 1/Q 2 read 

W 2 — M 2 ' 



A 



-2£n+[M'Z + 



n + A±, 



p = (l + 0n + 
p' + k = (l-i)n + 



M 2 . . W 2 — M 2 
(1 - O 

M 2 . W 2 -M 2 



n + 



2 ' 
A, 



Note that these expressions lead immediately to the relations 

2^ 



x B 



(12) 
(13) 
(14) 



(15) 
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and 

Ai = -(l-a(|A 2 |-|A 2 | mi „), (16) 
with the minimal value of the momentum transfer given by 

|A*U = + + W> - M 2 } = '° tf+ | "" t "- tf '' (17) 

1 — ? 1-xb 

Setting k — and W = M everywhere, one recovers the well-known formulas of the familiar 
case without pion. 

Finally, in order to quantify how the final nucleon and pion share the momentum with 
respect to the lightcone direction n, we further introduce a variable a such that 

k-n = a(l-0, p'-n=(l-a)(l-0- (18) 



3 Formulation of pion-nucleon parton distributions 
3.1 General parametrization 

In the case of ordinary DVCS or hard meson production, the factorization of the amplitude 
leads to nucleon matrix elements of lightcone operators. These nonperturbative objects 
are parametrized in terms of generalized parton distributions. For hard meson production, 
additionally the distribution amplitude of the final meson appears. 

In the present situation, where we take into account an additional pion, we arrive at 
matrix elements of lightcone operators with initial nucleon and final pion-nucleon state. In 
dealing with these objects, our first step is to provide a parametrization. Let us start with 
the isoscalar quark operator of vector type: 

^ e BA,p.n (Ar(p>a(A;) ^ ( _ An/2)#(An/2)|iv(p)) = J- U{p') F^VUip) (19) 

where we have introduced the Dirac matrices 

Ti=7 5 , r 2 = ^ 75 , r 3 = |-75, r 4 = i^ 75 (20) 

n ■ p M n ■ p 

with the definitions A = p' + k — p and p = {p' + k + p)/2. The constant qa is the 
isovector axial coupling constant, and f n is the pion decay constant normalized according 
to the experimental value of 93 MeV. Including these constants into the definition will be 
convenient in later calculations. Here and in the following, the insertion of the appropriate 
gauge links within the operators is always understood. The functions Hf \ that we shall 
call pion-nucleon parton distributions, depend on the momentum fraction x as well as on 
five further quantities that can be built from the vectors n, p, p', and k which we have at 
our disposal. We choose 

= H?\x,^A 2 ,a,t,W 2 ), (21) 



■5 



where in the present context of parametrization the definitions of the variables are 



In a similar way, we introduce the isovector even and odd distributions and 

— e 6A ^ n (iV7r a |V;(-An/2)^rV(An/2)|^) 

27T 

= ^j- U'Y^H^ + Qe^Ht^U, (23) 
■' 7r i=i 

which depend on the same set of arguments as H^°\ of course. 

The matrix elements of the quark operators of axial type differ from the previous ones 
only through the insertion of a matrix 75. Therefore, we can define in an analogous manner 
isoscalar distributions, 

J ^e fiA ^"(iW|^(-An/2)r/ 75 ^(An/2)|iV) = ^ E U'T^H^U (24) 

and isovector distributions, 

— e BXxp - n (N'Ti a \iP(-\n/2)^ 5 T b iP(\n/2)\N) 
2n 

= MfY^ U'Ta 5 (S ab Hl +) + te^H^U. (25) 
^ i=i 

Furthermore, the gluon matrix element can be parametrized as 

J ^e flAsp - n <iV'7r |F^(-An/2)F/(An/2)7i /1 7i I ,|JV) = E U'V^xH^U. (26) 

i=i 

At this point we have to comment on a corresponding parametrization that was previ- 
ously given by Bliimlein et. al. The authors came to the conclusion that five functions 
are necessary for a complete description. However, as demonstrated explicitely in appendix 
E of reference (TH], one can show that one of their functions can actually be reexpressed 
in terms of the others, i.e. the structures given in their parametrization are not linearly 
independent if one takes into account the Dirac equation. 

3.2 Reduction at the pion threshold 

Within the soft pion kinematics described in sectional all variables are required to be close 
to their values at the pion production threshold. Therefore, the exact threshold kinematics 
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can serve as a reasonable approximation e.g. for the calculation of scattering amplitudes. 
For this purpose, let us consider the threshold case and provide some useful relations. 
First, we recall that the invariant mass of the final nucleon-pion system is equal to 



th 



(M + m^) 2 . 



This is equivalent to the statement that the four-momenta k and p' are proportional, 

k = — p , 
M 



(27) 



(28) 



hence the number of independent invariants is reduced. In particular, we find that A 2 can 
be expressed through t and that a becomes a constant: 



M + m 

^th _ Jl/T 1 + m 7T5 



M 



«th 



M + m n 



(29) 



As a further consequence of the relation (|28l) . the Dirac structures involving e.g. T 1 = 
75 and T3 = ^75/M are no longer linearly independent. Therefore, the number of ttN 
distributions can be reduced by two. For example, in the case of the vector distributions 
Hi, we arrive at 



^e sx ^ n (N'n a \^-Xn/2)^(\n/2)\N) = ^flJ' 
2-K Mf n 



rr(0) , M i rr(0) 
lth + 2th 



l 5 r a U (30) 



and 



dA 
— ( 

2vr 



fi Xx P- n (N'n a \ij(-Xn/2)^T b i;(\n/2)\N) 



75^, (31) 



with the threshold pion-nucleon distributions H lth and H 2t h given by 



if. 



lth 
(0,±), 



'2th 

50U 

gluon operator 



f^fo e, A t 2 h , « th , t, ^ 2 h ) + ^iff £, A t 2 h , a th , t, , 
Hi°' ±] (x, C A t 2 h , a th , t, ) + ^H^ix, £, A t 2 h , a th) t, Wf h ) 



(32) 
(33) 



Analogously, we define #1°^ and for the axial operators and and -ff^th f° r the 



3.3 Moments of the ttN distributions 

The moments of ordinary nucleon GPDs are polynomials in the skewedness variable, where 
the coefficients give the nucleon form factors of the corresponding local twist-2 operators 
[El EI]- For example, taking the first moments of the quark GPDs yields 

J dxH(x,£,A 2 ) = F 1 (A 2 ), J dxE(x,£,A 2 ) = F 2 (A 2 ), (34) 
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dx H(x, £, A 2 ) = Ga(A 2 ), 



dxE(x,£,A 2 ) =G P (A 2 



(35) 



where Fi and F 2 are the Dirac and Pauli form factor and Ga and Gp the axial and 
pseudoscalar form factor, respectively. Concerning the second moment, we have e.g. 



Jdxx # (5) (x,£,A 2 ) + i# (G) (x,£,A 2 



dxx 



^)(x,e,A 2 ) + ^ G )(x,e,A 2 ) 



A(A 2 ) + C(A 2 ) (20 2 , (36) 
B(A 2 ) - C{A 2 ) (20 2 , (37) 



where A, B, and C are the form factors of the energy-momentum tensor. Similar poly- 
nomiality conditions hold for the moments of the irN distributions, their moments are 
polynomials in the variables £ and a. We shall demonstrate this now explicitely for two 
examples. 



3.3.1 First moment: The local limit 



The first moments of the irN distributions H^ '^ are related to the form factors of the 
matrix element which describes pion emission from the nucleon induced by the the local 
vector current (the hadronic ingredient of the pion electroproduction amplitude). For 
our purposes, the following parametrization in terms of pion emission form factors A4 is 
convenient: 



(N(p')ir a (k)\il)'f 
%4 



1 



i=i 



rM| 0) 

§ab A (+) + fi £ abc T c A (~) 



where the set of Dirac matrices is chosen as 

{rf, . . . , r£} = A", f, Y, 



(38) 



(39) 



(For a traditional parametrization we refer to Amaldi et. al. j!2j.) The form factors are 
functions of three independent invariants, e.g. A 2 , W 2 , and t. Current conservation reduces 
the number of independent form factors to six: 

(W 2 - M 2 )A 1 + 2A 2 A 2 + (W 2 + u-2M 2 )A 3 + 4MA i + 2(W 2 - M 2 )A S = (40) 



2A A + (W 2 - M 2 )A 5 + 2A 2 A 6 + 2{W 2 + u- 2M 2 )A 1 



0. (41) 



From the contraction of the matrix element (|38|) with the lightcone vector n, it follows 
that the first moments of the nN distributions are polynomials in A ■ n/p ■ n = — 2£ and 
k ■ n/p ■ n = a(l — £) = a: 



dx Hi = At- 2£A 2 + aA 3 



M / dxHo = A 



4- 



(42) 



-1 



-1 
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1 1 

fdxH 3 = A Fj - 2£A 6 + aA 7 , fdxH i = A 8 . (43) 



-l 



Note that the current conservation relations (|4T)|) and (|4*T| impose nontrivial conditions on 
the moments of the nN distributions. 

3.3.2 Second moment: The energy-momentum tensor 

The second moment of the 7rN distributions and is related to the form factors 
of the amplitude for pion emission induced by the energy-momentum tensor, which reads 

T» v = -ip^{D - Dy^ip + y —F pa F pa + F w F p u , (44) 

where the curly brackets denote symmetrization of the indices and D the covariant deriva- 
tive. We parametrize this amplitude as follows: 

„ 20 

(N(p') ix\k)\T^m P )) = j§fY. W) (45) 

^ i=l 

where the Dirac matrices are 

{If, . . . , T^} = {g^, ff, A»A", k»k u , pt"A">, p^k" } , A { »k u \ 

g^k 1 , p*f% A"A"A<, k"k u ^, p^A^k 1 , p^k v} k\ A^k v H, 
7 {^} ;7 {m A^, 7 ^F } , 

^7 {m , ^p v \ ^A v \ ^k u} } l5 . (46) 

As in the case of the vector current, the form factors Bi are functions of e.g. A 2 , W 2 , and 
t. From energy-momentum conservation, we have derived the following set of constraints: 

AB X + 4A 2 53 + (W 2 - M 2 )B 5 

+ (W 2 +u- 2M 2 )B 7 + AMB 16 + 2(W 2 - M 2 )B l9 = 0, (47) 



2(W 2 - M 2 )B 2 + 2A 2 E 5 + (W 2 + u - 2M 2 )B 6 + AMB lb + 2(W 2 - M 2 )B 18 = 0, (48) 

4(W 2 +u- 2M 2 )B A + (W 2 - M 2 )B 6 + 2A 2 B 7 + AMB l7 + 2{W 2 - M 2 )B 20 = 0, (49) 
AB 8 + 4A 2 B 10 + (W 2 - M 2 )B l2 + (W 2 + u - 2M 2 )B 1A + 2B W = 0, (50) 
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2(W 2 - M 2 )B 9 + 2A 2 B l2 + {W 2 + u - 2M 2 )B lz + 2B 15 = 0, (51) 
2{W 2 +u- 2M 2 )B U + (W 2 - M 2 )B 13 + 2A 2 B U + 2B 17 = 0, (52) 
(W 2 - M 2 )B 15 + 2A 2 B m + (W 2 + u - 2M 2 )B l7 = 0, (53) 

(W 2 - M 2 )B l8 + 2A 2 B l9 + (W 2 + u - 2M 2 )B 20 = 0, (54) 

which reduces the number of independent functions B { to twelve. The polynomiality con- 
ditions read 

i 

J dxx (Vj 0) + \Hf ] ^j =B 2 + 5 3 (20 2 + B A a 2 + £ 5 (-2£) + B 6 a + B 7 (-2£a), (55) 

-l 

i 

M j dxx (^ 0) + \h¥ ] ) = B 15 + 5 16 (-20 + B 17 a, (56) 



J dzx(#f + ^f } J = B 9 + B 10 (2£) 2 +B 11 a 2 + B 12 (-2£) + B 13 a+B u (-2£a), (57) 

-l 

i 

Jdxx^H^ + ^Hf^ = B 18 + B 19 (-2O + B 20 a. (58) 

-l 

Thus we see that from current conservation and polynomiality it possible to uniquely 
determine the form factors B^ from the second moments of the irN distributions + 
H\ G) /2. 

4 Soft-pion theorems for pion emission from the nucleoli 
induced by lightcone operators 

In the last section, we have presented a parametrization of matrix elements for pion emis- 
sion from the nucleon which is induced by twist-2 quark or gluon lightcone operators. Now, 
we turn to the calculation of these objects in the soft-pion region. For this, we rely on 
polology arguments, PCAC, and current algebra. The basic ideas are similar to those in 
the work of Guichon et. al. [BJ. However, while Guichon et. al. excluded small momen- 
tum transfer, we shall consider this region in our derivation as well. Actually, for certain 
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operators it was already shown in reference |13| that at small momentum transfer, the 
soft-pion theorems match the tree level results of a chiral perturbation theory treatment. 
This refutes opposite claim of ref. [14]. The results of (EJ and (TH] were confirmed also in 

ref. ma. 



4.1 Simple pion emission from the nucleon 

First, in order to introduce some notations and to demonstrate our approach on a very 
simple example, let us consider the emission of a pion from the nucleon, i.e. the amplitude 
M(N'ir\N) defined through 

M(N(p')7r a (k)\N(p))= lim *^^<j\^(p')|^|j\^(p)), (59) 

where $ a is the interpolating pion field 

d-A a 



A"; the axial current, 



$ a = V^r, (0|$'V) = 5 a \ (60) 

K = ^7,75 (61) 
and (only within this subsection) k denotes the nucleon momentum difference, 

k» = (p- p 'y. (62) 
If we relax somewhat about the on-shell requirement for the pion, we can also write 

(N(p)7i a (k)\N(p)) = (2n) 4 5(p' + k-p)M(N(p')7i a (k)\N(p)). (63) 

To derive the soft-pion theorem for the amplitude Ai(N f 7i\N), we start from the nucleon 
matrix element 

I u = ±r{N(p>)\Al\N(p)). (64) 

Jtt 

Applying the definition of the interpolating pion field (jHOI) yields the identity 

k v I v = timl(N{p')\<5> a \N{p)). (65) 
This equation is now investigated in the region where k is small, 

k ~ rriir ~ e 

i.e. we search for the leading contributions in e to both sides of equation (Jfiojl . 

The nucleon matrix element on the left hand side is parametrized in terms of the axial 
form factor Ga and the pseudoscalar form factor Gp, 



I l/ = ±(N'\A a u \N) = ±U' 

Jtt Jtt 



G A (k 2 ) lu -G P (k 2 )- kl 



2M 
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75yt/. (66) 



For small A; 2 , Ga approaches the axial coupling constant and Gp is dominated by the 
corresponding pion pole contribution, 

G A (k 2 ) =g A + 0(e 2 ), G P (k 2 ) = + 0(l/e). (67) 

K — m n 

The insertion of these relations and the use of the Dirac equation yield 



,2 



WI V = -jJ^U'^-^U + 0{e 2 ). (68) 

The right hand side of equation (|HH|) has an explicit factor of m\ ~ e 2 . This small 
factor can be compensated only by the pion pole contribution, 

fim 2 (N(p')\$ a \N(p)) = Qm 2 6 - M{N{p')ir a (k)\N{p)) + 0(s 2 ). (69) 

k — 

So finally, the equation (j65| leads to 

M(N'7r a \N) = U'^l/hX-U + 0(e 2 ), (70) 

Jit ^ 

from which in turn the Goldberger-Treiman relation emerges, if we parametrize AA(N'tt\N) 
in terms of the pion-nucleon coupling constant g n NN as usual. The line of argument that 
lead to this well-known result is now generalized to the case when a lightcone operator is 
present. 

4.2 Pion emission induced by lightcone operators 
4.2.1 Isovector quark operator of vector type 

Now we turn to the pion emission induced by lightcone operators. First, we deal with the 
isovector quark operator of vector type which we shall refer to as O b (X) in the following: 

O b (X) = $(-An/2) yirV(An/2). (71) 

We recall that in terms of local operators it reads 



K ' ^ ml V 2 

■m— n \ 



^(0) [n-(d - d)]"YrV(0). (72) 



- {-) <^"i'"^ - or- 1 -- 

The derivation of the soft-pion formula for (iV'7r|0 6 (A)|./V) starts from the following 
object denoted I u : 

h = jj ^ze^{N{ P ')\T[Al{z)0\X)]\N{ P )), (73) 
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where T stands for the time-ordering prescription 

T[A a u (z) 0\\)] = 9(z ) A a v {z) 0\\) + 9(-z ) O b (X) A a u (z) (74) 
Let us investigate the behavior of the matrix element l v in the soft-pion region 

k ~ ~ 6. (75) 
In this region, several momenta are close to the pion or nucleon mass shell: 

k 2 = ml + 0{e 2 ), W 2 = [p' + kf = M 2 + 0{e), u = {p-k) 2 = M 2 + 0{e). (76) 
Generally, for a Green function of n Operators, 

G x ... n {(3\a) = f d A Zl . . .d^e^-^"-^- 1 ^^ 

(77) 

we expect the following behavior when \a) = \(3) = |0) and the momentum p = Pi + . . - +p r 
reaches the mass shell of a particle with mass m and quantum numbers s: 

Gl„»(0|0) ^ + - + ^ w V^.j^C+ufe^O) (78) 

s 

(see, e.g., Chapter 10.2 in the textbook of Weinberg |IH]). The matrix element l v can be 
related to a Green function of the type Gi... 4 (0|0) according to the usual LSZ formalism. 
However, in our case we deal with the additional problem that several momentum com- 
binations become nearly on-shell at once, as shown in the equations (f76l) . We proceed by 
adding up these different pole contributions, each of which is of the type (|T%|) . and subtract 
appropriate terms where necessary to avoid a double counting. In this way we obtain the 
preliminary statement 

F = Kik) + I U N {P ' + k) + J V*) + °( £ °)> (preliminary) (79) 
where I\ h -, denotes the pion pole in the variable k 2 , 

F ^ Uo\A» a (0)\n c (k)}— ^(N(p')n c (k)\O b (X)\N(p)} = F <k) = 0(l/s), (80) 
Jtt k — m„ 

^N( P '+k) ^ ne nu cleon pole in the variable W 2 with pion-pole subtraction, 

w 2 ^m 2 1 



TV TV 



T l 

1 TT(k) r 
Jtt 



(N(p')\A» a \N a (p' + k)} 
-(O|^> c (*0) , 2 6 j M(N{p')7r c (k)\N a (p' + k)) 



(81) 



X W 2 -Mi {N ° iP ' + k )\° b W\ N (P)) = J V+*) = °W e )> ( 82 ) 
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and the nucleoli pole in u with pion-pole subtraction, 

[ {N a (p-k)\A™(o)\N( P )) 



l v -K(k) ^ UN(p')\O b (\)\N a (p-k))- l? ' 



(0|^ a |7r c (fc))— ^^(iV^p-fc) 7r c (A;)|iV(p)) 



(83) 

'V*) = °(Ve). (84) 



In the previous formulas, c denotes the pion's isospin and a a combined nucleon spin-isospin 
index. Summation over quantum numbers of intermediate states, such as c and a here, are 
henceforth understood. 

An additional pole that we have not mentioned yet has to be considered in the kine- 
matics of small momentum transfer, t ~ e 2 , since then 

t = ml + 0(e 2 ). (85) 

It has the form 

r d 4 z e Bk - z (0\T[A va (z) O b ] \ir d (p - p')) 
B 



TV TV L " l ir J - 

1 ~ 1 n(k) > ~T~ 



iO\A^\^k)) w ^(7r%k)\O b \Ap-p')) 



(86) 



K M{N'-K d {p- V ')\N)=I^_ pl) , (87) 



t — mi 



where a 7r(fc)-pole subtraction has been performed. Therefore, to cover the whole region 
of momentum transfer —t < M 2 , we modify equation (fT^l) by adding this contribution, 

I" = K(k) + ^N(p'+k) + I U N{p-k) + K(j>-p') + 0( £ °)i (88) 
always keeping in mind that I^ p _ p ,^ is suppressed automatically as —t grows: 

r _ f 0(l/s) t ~ e 2 
Up-p') - \ O(s ) -t > e 2 ■ 

This closes the discussion about the pole structure of the matrix element I v . 

Next we turn to a simple identity for I u that follows when we use the definition of the 
interpolating pion field (|60|) : 

k u I u = Q j d 4 ze m - z S^^(N'\[A%(z),O b (\))\N) + (iV'|T[$ a (z) O b (A)] |iV) | . (90) 

Let us exploit this identity in the soft-pion region (f75|) . According to our previous demon- 
strations, the leading contributions to the left hand side are simply obtained by contracting 
k with the poles in (jHE} - The contraction with the Tc(k) pole yields 

kvK [k) = -^-^(N'7r a (k)\O b (X)\N) + 0(e), (91) 

7T 
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where we have used the definition of the pion decay constant, 

(0\A a u \n\k))=8Uk v 5 ab . (92) 
For the sum of the nucleon poles we obtain 

+UT^(p',p- k, A) r b ^^lj-h^U + 0(e), (93) 

where we have applied 
k v 



f (N(p 2 )\A:\N{p 1 + k)) = --^_U(p 2 )^ l5 T a U{p 1 ) + 0{e 2 ) 

Jit K — m 



and 



9A 



M(N(p 2 )n a (k)\N( Pl )) = U(p 2 )^ l5 r a U( Pl ) + 0(e 2 ), 



(94) 



(95) 



as discussed in the previous subsection. Moreover, we have introduced the parametrization 
(N(p 2 )\O b (\)\N(p 1 ))=U(p 2 )r( v \p 2 ,p 1 ,\)T b U(p 1 ), (96) 
where in terms of generalized parton distributions one usually defines 



r (y) (P2,Pi,A) = / dx e~ SXxiP2+Pl) - n 



i/2 



H (VU | E ( V) Mn,P2-Pi) 
1 2M 



(97) 



with 

and analogously for E^ v \ Further, we mention that it is useful to rewrite the contraction 
with the 7i(p — p') pole in the form 



k T v 

IS 7T 



(p-p') 
6 

u 

k 



j- \ d 4 ^e 6fc - 2 { ( 5(^)(0|KW,Ol|vr d (p-p0)+^(0|T[^^)O 6 ]|^(p-p0)} 



6 



^(0\A" a \lT C (k)) 19 - A 7C C (k)\O b \7T d (p-p')) 



fn 



k 2 — m 2 



B 



t — ml 



M(N'7i d ( P -p')\N) (99) 



l.(0\[Q^),O b ]\n d (p-p')) + Qml J d 4 ze &k - z (0\T[^(z)O b }\n d (p-p')) 



k 2 — m 2 



(n a (k)\O b \n d (p-p')) 



l! U'j^(p , - P "h 5 T d U + 0(e), 
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(100) 



where Q§ is the axial charge, 

Q a 5 (z )= [ d 3 zA a (z). (101) 



So much for the left hand side of the identity ([HOI) . Now, what are the main terms on the 
right hand side? 

First, we turn to the commutator term. It can be approximated by neglecting the soft 
momentum k, 

6 f d^-'^(N , \[j^(z) 1 O b (X)]\N) = ^(N'\[Qt(0),O b (X)]\N) +0(e). (102) 

For the calculation of the commutator [Qf(0), O b (A)], we insert O b (X) in the form (f72|) . 
Then we have to deal with objects of the type [<2f(0), (n-d) m ip(0)} and [Qf (0), (n-<9) m ^(0)]. 
For m = 0, these commutators are well-known from the transformation of the fields under 
the axial part of the chiral rotation: 

[Qt(0), m\ = -75y^(0), [O|(0)^(0)] = -^(0)75y- (103) 
In the general case m > 0, we can rewrite e.g. 

[Ql (n ■ d) m i,} = jr ( m k \ {-\)\n ■ dr~ k i(n d ) k Ql # (104) 

The derivative of the axial charge is given as 

doQt(zo) = j d 3 z^z) l5 {r a ,m}^z), (105) 

where rh = diag(m u , ma) is the quark mass matrix. According to the Gell-Mann-Oakes- 
Renner relation we have 

m u ,m d oc ml ~ e 2 , 
hence we can neglect the terms that involve derivatives of Q^: 

[Ql (n ■ «9)"V] = (n • d) m [Ql ^ + 0(s 2 ) (106) 

and similar for [Q 1 ^, (n • d) m ijj}. In this way we obtain finally 

[Qjf(0), 0\\)] = &e abc O c 5 (\) + 0{e 2 ) (107) 

with 

O c 5 {\) = ip{-Xn/2) ifi lb T c 4){\n/2). (108) 

Next, we discuss the second term on the right hand side of |90|) . It is accompanied by a 
small factor m 2 ~ e 2 , so again we focus on the pole contributions which are of order 1/e 2 . 
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Figure 2: Illustration of the soft-pion theorem for the matrix element 

(N(p')7T a (k)\O b (X)\N(p)), where O b (X) is the lightcone operator ^(-An/2) f£rV(An/2). 
The blobs denote insertion of the indicated operator while the pointlike vertices represent 
the standard pseudovector pion-nucleon coupling. 



These are the pion pole in k 2 and, for small momentum transfer, the pion pole in t. Again 
taking into account a double-pole subtraction, we thus get 



8ml J d 4 ze Sk - z {N'\T[$ a (z) O b ]\N) 



, , ^(N'-K a (k)\O h \N) + 6m' / d i ze Bkz (0\T\^ a {z)O b }\7c d {p-p')) 
k — m w 



+ 



m 



*—(7t\k)\0\\)W\p-p')) 



6 



,,_ mV ^ — r r vj , _ /;4 ^(iW^-p')|iV) + 0( £ ). (109) 

Now we can collect all approximations and insert them into (f90l) (note that the terms 
which involve the time-orderd product T[§ a (z), O b ] cancel). Then we can solve the equation 
for the matrix element (N'7r a \O b \N) and arrive at the soft-pion theorem 

(N'7r a (k)\O b (X)\N) 



-abc 



(N'\O c 5 (\)\N)-(0\OI(\)\n d (p-p>)) 



t ~ ml 2U 



- p"Hr d U 



2/. 



+ (7r a (k)\O b (X)\n d (p - jfy-^U'Jg-tf - f) l5 r d U + 0(e), 

see also figure [2j The pion matrix elements in this formula are parametrized as 

i 

(0\O c 5 (\)\7i d {p-p')) = 2$U(p-p')-n5 cd J dxe- eXx{p - p ' yn/2 Mx), 

-i 

where <j) n is the pion distribution amplitude, and 

(n a (k)\O b (\)\7t d (p-p')) 

l 

(p — p' — k) ■ n 



U 
(110) 



Qe abd (k + p-p')-n dx e -ZM*+p-p')-n/2 R \v) 



x. 



-i 



(p — p' + k) ■ n ' 



+0(8*), 
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In the last line, we have neglected the momentum transfer in the argument of the pion GPD 
because the n(p — p') pole terms only contribute in the region where t ~ A 2 ~ e 2 . 
For a parametrization of the nucleon matrix element of 0|(A) in terms of GPDs see below, 
formulas ljT33j) and (11531) . 

For reasons of consistency, it is advisable to check that the n{p—p') pole contributions in 
the soft-pion formula is really negligible for — t ^> e 2 , as we have promised at the beginning 
of the calculation. In this region of moderate momentum transfer, we have to distinguish 
two cases which are both kinematically allowed: 

(p — p') ■ n ~ k • n ~ e or {p — p') ■ n ^> k ■ n ~ e. (HI) 

In the first case, both pion pole terms are individually suppressed because the pion matrix 
elements are proportional to combinations of (p — p') • n and k ■ n. In the latter case, the 
combination of both terms is small because 

1 

(k + p - p) ■ n f dx e -&Mk+P-P>>n/2 H ( V ) ( (p-^-fc)-n Q \ 
J V {p-p' + k)-n J 

-l 

i 

(p-pO-n^-n-e (p _ pl) . n J dx e ~*>*to^l* H V)fr 1, ) + 0(e) (112) 

-1 

and a soft-pion theorem for the pion GPD [12] reads 

HS n (x,l,0) = 2<f> v (x). (113) 

So we have confirmed that in the region of moderate momentum transfer the ir(p — p') 
pole is negligible. Therefore, our result coincides with that one of Guichon et. al. [H], who 
restricted themselves to the region — t 3> e 2 . 



4.2.2 Isocalar quark operator of vector type and gluon operator 

Soft-pion theorems similar to the previous one can also be derived for the isoscalar quark 
operator and the gluon operator. For the isoscalar quark operator, we get 



(N'n a (k) \if>(-\n/2)fol>(\n/2)\N) 



^ J ^T^\p> + k,p, A) - W\p\p - k, X)^fHr a 



u 



+ (n a (k)\4,(-\n/2)^(\n/2)\n c (p - p'))-^U'^ - f)j 5 r c U + 0(e), (114) 

where can be obtained from the definition of in equation (|9T|) by replacing the 
isovector GPDs and E^ v > with the isoscalar ones, H^ s > and E^ s \ For the gluon 
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operator, we obtain 



(N'it a {k)\n^{-\n/2) F pi/ (Xn/2) n u \N) 



Hr J ^±^T^( P ' + k,p, A) - T^(p',p - k, A)£±4W 



2p> ■ k 



2p ■ k 



U 



+ (it a (k)\n,F"\-\n/2) F pv {\n/2) n»\n c (p - p >))-JL_u> l^ty _ f) l5 r c U 



+0(e), 



(115) 



where 



(116) 



with arguments of the gluon GPDs and E^ G ' as in (|HH1) . Because of the isoscalar 

nature of the two operators considered here, commutator terms do not appear in the 
results. Moreover, we remark that the discussion of the pion pole terms at moderate 
momentum transfer which we presented in the isovector case, can be repeated here in a 
similar way using the following soft-pion theorems for the gluon and the isoscalar pion 
GPD [Ilj: 

H {s \x, 1,0) = 0, H (G \x, 1,0) = 0. (117) 
So again, for —t ^> e 2 the pion pole contributions to the soft-pion theorem are negligible. 

4.2.3 Isovector operator of axial vector type 

Next we turn to the derivation of the soft-pion theorem for the matrix element of the 
operator 

0\{\) = TP{-\n/2)^ B r b iP{\n/2). (118) 

For this purpose, we define 



d *z e ^(N(p')\T[Al(z)O b 5 (\)]\N(p)) 



(119) 



and obtain the identity 



k u I v = H J d 4 ze 6fc - 2 ^^(N'\[A a Q (z),O b 5 (\)]\N) + m^T^) 5 b (A)]| AO } • (120) 

Again, we make soft-pion expansions of both sides when k ~ ~ e. At this, we have to 
consider pion poles in k 2 and A 2 = [p' + k — p) 2 and nucleon poles: 



I" — K{k) + K(A) + ^N{p'+k) + -^TV(p-fc) + 0( £ °)- 



(121) 
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The ir(k) pole is given as 



|(0|^(0)|tt c (A;))-^— (N'^(k)\O b 5 (\)\N) 

Jit K — III 



-k V 



k 2 — ml 



(N'n a (k)\O b (X)\N) = r <k) 



(122) 
(123) 



For the 7r(A) pole we get 

1 ~ 1 7T(k) > - 



-^<iV'|^ a (0)|iV7r c (A)> + - - M(N'Tr a (k)\Nir c (A)) 

Jit k — m w 



6 



x- 



(7r c (A)|O b (A)|0) = /: (A) , 



(124) 



A 2 — m 2 

where M(N'7r d (k)\N7r c (A)) denotes the pion-nucleon scattering amplitude, 

(iVV(A;)|iV7r c (A)) = (2tt) 4 5(p' + k - p - A) M(N'ir d {k)\Nir c (A)). (125 
This pion pole is non-negliglible only if —A 2 < m 2 ~ e 2 . The nucleon pole terms are 



JU JU JU 

1 ~ I n(k) ~ 1 tt(A) 



±(N'\A™\N a (p' + k)) + — - 2 M{N'ix a {k)\N a {p' + k)) 
Jit k — m n 



B 



W 2 — M 2 



x 



{N a (p' + k)\O b 5 \N) - M(N a (p' + k)\Nn c (A)) 



B 



A 2 



mi 



■(n c (A)\O b 5 \0) 



2/ w 2p' • fc 



rr) ( y + ^,A)^ + ^A 7 ^^ A )'^l > 



2/. 



A 2 
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and 



p _ p p 



(N'\O b 5 \N ff {p - fc)) - M(iV'|iV a (p - A;)tt c (A)) 



6 



A 2 — m 2 



(^(A)\O b 5 \0) 



6 



x 



u — M 2 



^(N a (p - k)\A» a \N) + — - 2 M{N a {p - k)n a (k)\N) 
Jtt k — m w 



U' 



-2p ■ k 2f w 



+0(e°) = r N(p _ k) , 



where we make use of the parametrization 

(N(p 2 )\O b (\)\N( Pl )) = U(p 2 )T { P(p 2 ,p 1 ,X)T b U(p l 



(126) 



(127) 



(128) 



(129) 



(130) 



(131) 



(132) 
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with 



i 



1 (2M) 2 



75- 



(133) 



Note that in contrast to the previous case of the vector type operator O b (A), a 7r(A) pole 
subtraction has been necessary to obtain the correct nucleon pole terms ()128ft and (|131|) 
without double counting. 

From the soft-pion expansion for I u that we have determined now, we immediately 
obtain the left hand side of the identity (|12()jl . k v l v . The right hand side gives 



8 j d 4 ze 

-abc 



fife- 



( S(Zp) 



(]V'|[^(z),0 5 b (A)]|iV)+m2(iV'|r[$ a (z)0 5 b (A)]|iV) 



£ 



■(N'\O c {X)\N) 



mt 



+ 



(N'\A a u \Nn c (A)) + 



k 2 — m 2 
m : 



k 2 — ml 



(N'n a (k)\O b 5 (X)\N) 
■M(N'n a (k)\Nn c (A)) 



fi 



A 2 — ml 



(n%A)\O b 5 (X)\0) 



+0(e). (134) 
In this way, we obtain from (|120p the following first version of the soft-pion theorem: 



(N'n a (k)\O b 5 (X)\N) 

g-abc 

■(N'\O c (X)\N) + M(N'TT a (k)\Nir c (A)) 



fi 



fn 



A 2 



mt 



>y{A)\o\{xm 



2U 



2p' ■ k 



>(p +k,p,X)r + — A75T 



-U' 



T^(p',p-k,X)r b + ^ l5 r' 



2U 



% 

2/tt 

<tt c (A)|O*(A)|0> 



<tt c (A)|O*(A)|0> 



A 2 — ml 



U 



A 2 



mi 



j> + M^g A 
2p ■ k 2f n 



+0(s). 



(135) 



To further simplify this expression, we need the pion-nucleon scattering amplitude 
which appears in the vr(A) pole terms in the region of small momentum transfer. We can 
derive the soft-pion theorem for this amplitude from the identity 

A M 



(N'7r a (k)\A c J0)\N) = -fim 2 (iVV(£;)|$ c (0)|iV). 



(136) 



For the matrix element (N r ir a \A c \N) on the left hand side, we can use the soft-pion theorem 
(|135j) since 

^(0) = 20 5 C (A = 0). (137) 
The right hand side of (|136|) is dominated by the pion pole in A 2 , 

-,2 



-fim:(iVV a (£;)|$ c (0)|iV) = M(N'n a {k)\N7T c {A)) 



mz 



A 2 — ml 



0(e 2 ). 



(138) 
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m = v 




Figure 3: Leading contributions to the pion-nucleon scattering amplitude 

M.(N{p') 7r(A)|iV(p) 7r(k)) for soft pion momenta A and k. The pointlike vertex with 
four particles attached corresponds to the Weinberg coupling. 



Thus, we can solve equation ()136|) for the pion scattering amplitude and obtain the well- 
known result 

^acd 

M(N'<K a {k)\N<K c (A)) = -—U'{ft + fi)r d U 

see also Fig. El If we insert this expression into equation (11351) . we find that the double 
poles are canceled, and the final version of the soft-pion theorem reads 

(N'n a (k)\Ol(X)\N) 

■ -UT^(p',p, A) r c U - £ —U\$ + fi)r d U- -<tt c (A)|O*(A)|0> 



U ' 4/2 " ^ A* -ml 

+0(e), (140) 
see also figure EJ 



Moderate momentum transfer In the region of moderate momentum transfer, i.e. —A 2 
^> e 2 , the pion pole term in the soft-pion theorem (|140|) is negligible as it should since ac- 
cording to the Dirac equation for the nucleon spinors we have the suppression factor 

U'(jk + fi)U = 2U'jkU = 0{e), (141) 

while the denominator A 2 — m 2 is no longer small. Therefore, in this region the soft-pion 
theorem again agrees with the corresponding one in Guichon et al. [HJ. 
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£ abc O c /f n 



X 



o 



+ 




o 



+ 




A 



Figure 4: The soft-pion theorem for the matrix element (N(p') 7c a (k)\Ol(\)\N(p)) , where 
O b B (X) is the operator ^(-An/2) ^T b ^{\n/2). 



4.2.4 Isoscalar operator of axial vector type 

In the case of the isoscalar axial quark operator, matters are simpler because commutator 
and pion pole terms vanish and we are simply left with the nucleon pole terms: 



(N'ir a {k) \ij(-\n/2)^ 5 i;(Xn/2)\N) 
ft + M 
2p' ■ k 



where 



2U 
+0(e), 



r ( 5 s) (p 2 , Pl ,X) = I dxe- SXxn - (pi+P2)/2 



Ur J ^T^ S \p' + k,p, A) - rf V,p - k, A) 



]/> + M 
2p-k 



u. 



1 2M 



7s- 



:i42) 



(143) 



-i 



4.3 Application to the pion-nucleon distributions 



From the soft-pion theorems for the matrix elements of twist-2 lightcone operators, equa- 
tions (|110|) . ()114p . (|115j) . ()140|) . and ()142|) . we arrive at the corresponding results for the 
7iN distributions Hi and Hi by Fourier transformation with respect to the light-cone co- 
ordinate A and an appropriate decomposition of the Dirac matrix structure. These results 
are presented in the following at the pion threshold. The general (non-threshold) results 
are given in appendix iBl 

The soft-pion theorem ()110p leads to 

<W,*) = -^f^(^,AV^g^^(^,A 2 ) (144) 

= " g " /,( ^ 2 g; £ ;; ,( ^ A2)i ™ 
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and 

H&(x,t,t) = ^[EW(x,to,t)-EM(x,t ,t)} 
9a 

2M 2 -£ K ; 2(2M 2 -£) V ; 

+^~A V) (^-> o) ■ " - N) (146) 
m 2 - £ \p t -n pt-n J 

H (-) (rff) H(x,Zo,t) t[g^)(x,e,A 2 ) + ^)(x,e,A 2 )] 

^th^,*) " ^ 2(2M 2 -£) • 

Here, A 2 is considered as function oft according to the threshold requirement (|29j). The 
variable £o is defined as 

6 = fe4^. (148) 

its threshold value is related to £ via 

th (2M + m 7r )e + m 7r 
e ° " 2^1^ + ^(1 + ' ( 49) 

The average momentum in the £ channel is referred to as pt, 

p — p' + k , 
Pt = P — 2 • (150) 

At threshold, p t ■ n can be expressed in terms of £: 

th c m n (l -f) 
Pfn = i + —— . 151 

M + m n 

Further, E^ e denotes the pion-pole contribution to E^ and is given in terms of the pion 
distribution amplitude 6„ as [T8l H9l l20l ITT] 

SSZM,t) = ^^U(f)<H6,-M). (152) 

m n ~ 1 SO \S0/ 

Next, we turn to the 7tN distributions for the soft-pion theorems (|114|) and (|115|) . They 
read 

+^7 fff G) (— • — • o) »fe ' " " W) (153) 
m 2 - £ \pt-n pt-n J 

ff (o,c )(r , f) t[H^(x^A')+E^(x,^)] 

^ 2th ^,^«J - 2 (2M 2 - £) ' 1 j 
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where the first superscript refers to the isoscalar quark distributions and the second one 
to the gluon distributions, respectively. 

Finally, we give the results related to the soft-pion theorems (jl40jl and (j!42jl . again 
using a combined notation for the isoscalar (0) and isovector even (+) GPDs: 

>2\ 



and 



H 



(o,+), 

lth 



H. 



(o,+), 

2th 



2M 2 -t 1 ' C ' ' 2(2M 2 -t) 

AM 2 -t 



2(2M 2 - t) 



m s > v \x,£,A 2 



E^ M 2^(1-0^, ^ 



fl'A 2M 2 - t 

t£E( v \x,£, A 2 ) 2Mm 



2(2M 2 -t) 



C/.4 



(155) 
(156) 



-^(f) ^-N) ( 157 ) 



2(2M 2 - t) 

4M2 - J gw (j , t , A . ) _^°(».&.«)+g <v) (».fa«). (168) 



4.4 Results for the moments of ttN distributions 

In section I3.3| we have described how to obtain pion emission form factors of certain 
local operators from the moments of the ttN distributions according to the polynomiality 
property. Now that we have derived the soft-pion theorems for the nN distributions, we 
can easily read off these form factors after taking the corresponding moment. In the cases 
where the results are known, this procedure provides a check of the previous calculations. 



4.4.1 First moment of Hf^ 

The soft-pion theorems for pion emission induced by the local vector current, i.e. for the 
form factors Ai (i — 1, . . . , 8) that we have defined in formula (|38j) . are given in appendix 
ICl For a discussion, we restrict ourselves in the following to the pion threshold where the 
number of independent form factors is reduced to two. The conventional quantities for 
such a description are the transversal and longitudinal s wave multipoles E^ '^ and L( 0,± ). 
In the center-of-mass frame where p ' = k = and p = —A, these threshold multipoles can 
be defined using the spatial components of the matrix elements 



(N f (p',S')n a (k)\^A \ \i/>\Ni(P,S))=Z< 



T-iM(O) a 
1 S'S 'ft 



TPs 6 $h + T^-'Qe abc r 



S'S 



'fi 



namely 



6 s ' s 
2 s ' s 



E^ffps + (Lq+ - E { 



(0) 



0+, 



\p\' 



F {±)-* , (T (±) F (±h PP-<rs>s 
A)+ a s's + l^o+ - A)+ ) — j^p — 



(159) 
(160) 
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where e is the electromagnetic coupling constant and a = (a x ,a y ,cr z ) are the Pauli spin 
matrices. 

According to these definitions, we obtain for the relations between these multipoles and 
the 7iN distributions at threshold 



E, 



(0) 



^i?(A 2 ) 
2(2M 2 - 1) 

tGZ\A*) 



E, 



(+) 



-ec J dxH^(x,£,t) = ec 



-i 



0(e) 



0{e) 



E 



(-) 



-ec J dxH^(x,£,t) = ec 



2(2M 2 - 1) 



FP(t) + 



2(2M 2 - 1) 



0(e) 



(161) 
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with the kinematical prefactor 
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+ 0(e) 
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&a V4M 2 - 1 
I6nf n M + m-x 



The functions 

G M (A 2 ) = F 1 (A 2 ) + F 2 (A 2 



G E (A 2 ) = F 1 (A 2 ) 



A 2 



(167) 



(168) 



are the magnetic and electric form factors of the nucleon, with the superscripts S and V 



AM 2 ' 
: sup 

indicating the isovector and isoscalar combination, respectively. F^ = G^/9a denotes 



the isovector axial form factor normalized to unity. Figure El shows a plot of the multipoles 
E 0+ and L 0+ as function of the invariant momentum transfer A 2 . 

The multipoles E^ + ^ and turn out to be quite small, so that the corrections 

which are neglected in the soft-pion theorems become important for these quantities. In 
the conventional units that are also used in figure these corrections could be as large 
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Figure 5: The transverse and longitudinal s-wave multipoles as function of A 2 = (1 + 
m 7T /M)t + m 2 . The curves up to —A 2 = 1 GeV 2 correspond to the soft-pion theorems 
(ITFTT) to (IT661) . The curves up to -A 2 = 0.1 GeV 2 show the one-loop ChPT results by 
Bernard et. al. \'2'2\ as given in their formula (5.1). The references for the nucleon form 
factors which have been used for this plot can be found in appendix For the pion 
electromagnetic form factor we have used the monopole form F^A 2 ) = (1 — A 2 (r 2 )/6) _1 
with (r 2 ) « 0.45 fm 2 , see references [23l l24"j . 



as cm n /M ~ 3 x 10~ 3 /m 7r +. The situation is different for the multipoles Eq + and Lq + ; 
they reach sizable values in particular at low momentum transfer, and corrections become 
relatively small. 

Since the first use of soft-pion theorems in pion-electroproduction that can be dated 
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back to the work of Nambu and Shrauner in 1962 [25], a lot of progress has been made 
in calculating corrections to the early soft-pion results. Within the framework of chiral 
perturbation theory, pion-electroproduction has been computed including the one-loop 
level |22j. The resulting threshold multipoles are also shown in figure El Unfortunately, 
these achievements of ChPT are limited to the region of very small momentum transfer, 
—t <C M 2 . If we expand the soft-pion results ()16H) to (|166|) for such small —t, we get 
correspondence with the leading-order results of ChPT. 

4.4.2 First moment of 

Further, we would like to discuss the form factors of pion emission induced by the local 
isovector axial current. They can be obtained from the first moments of the rcN distri- 
butions h\ . Again, we restrict ourselves to the threshold. Then the matrix element is 
parametrized in terms of three s-wave multipoles Lq+ ■> Hq+ > an d M)+ : 

(N f (p', S>) n^lA^mp, S)) = \{5 ah 5 f ,{f^ ] + fie^T^), (169) 

where in the center-of-mass frame (p ' + k = 0) one defines for the time and space compo- 
nents of T M 

f°Jp = 16tt(M + 7^X45 + A #£Vs's, (170) 
T$> = 16tt(M + (-M+Ss's + M$S as '** P ) (171) 

(see [221 and original references therein). The relations between these multipoles and the 
moments of theshold tiN distributions are 

4? = c j Ax (M+p^ mxAt) + 2M ?^ m*,(,t)) (172) 

-1 

1 

= ^/^(-^(^l^)-^^^^^*)) (173) 
-l 

i 



-i 



where c is the same kinematical factor as in definition (|167|) . Note that in taking the 
moments of the ttN distributions at a particular skewedness £, in general one has to 
integrate first and then to insert the value of £ to avoid unphysical regions in the £-t 
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plane. After the insertion of the soft-pion theorems for the functions Hah, we obtain 

~Gff(t) {AM 2 -t)GP{t) 



M, 



(-) 



o+ 



M, 



(+) 



9a 
16irf n 

9a 



M 2 



9A 



2{2M 2 -t) 
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(175) 
(176) 

(177) 



while Lo+? an d #0+^ are vanishing within the accuracy of soft-pion theorems. At small 
momentum transfer, these results agree with the leading order expressions of the ChPT 
calculation in reference EI 



4.4.3 Second moment of Hf ] + H^> /2 

As described in section l3~3| the second moment of the irN distributions and give 
the form factors B± for soft pion emission induced by the energy-momentum tensor. The 
resulting soft-pion theorem is 
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l5 r a U + 0{e) 



(178) 



where A = A{A 2 ) and so on for the other nucleon form factors B and C. For the determi- 
nation of Bi and B 8 from the nN distribution moments, it has been necessary to use the 
current conservation relations (jITjl and (jHUjl . because these form factors appear together 
with the metric tensor g^ u which vanishes after contraction with the lightcone vectors n^n" ' . 
The remaining six current conservation constraints are fulfilled, so that 



A^N'tt^T^IN) = 0(e) 



(179) 



5 Hard pion production with additional soft pion 

As an application of the presented soft pion theorems for the nN distributions, we consider 
the process of hard 7r + production off the proton with soft pion emission. The two possible 
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reactions are 

ll{q) + p(p, S) -> n(p', 5') + 7r°(fc) + vr+(g'), (180) 

and 

ll{q) + p(p, 5) - p(p', 5') + Tr-(fc) + 7r+( g ')- (181) 

We shall calculate the longitudinal differential cross section as well as the transverse spin 
asymmetries for these processes, in comparison to the corresponding pure process, namely 

72(g) + p(p, S) -> n(p', 5') + 7r + (g'). (182) 
5.1 Longitudinal cross section 

First, we recall the definition of the differential longitudinal cross section of 7r + production 
without soft pion emission, reaction (|182j) : 



dVf = \^Sl^M±MM (183) 
2(s-M 2 ) v /A(s,M 2 ,-g 2 ) 

where d 2 $ is the two-particle phase space volume, 

d 2 <% + p; q>, p>) = ^l^l(2^S(q + p-q'-p'), (184) 
and A denotes the conventional kinematical function 



A(x, y,z) = a/(x - y - z) 2 - Ayz. (185) 

Introducing the angle 0, which is the azimuthal angle of q ' with respect to the direction 
of q in the center-of-mass frame, and t = (p — p') 2 , one obtains for the cross section 



|2 



2.» lA^Wl 

(186) 



dV[ nj _ \M\ 



d(j)dt 32n 2 (s- M 2 )A(s,M 2 ,-Q 2 ) 
The amplitude is given through the matrix element 

MP = e{n(p')n + (q')\e L .J\p(p)) (187) 
of the electromagnetic current 

J = (I + ^] 1>, (188) 



.6 2 t 

where El is the longitudinal polarization vector, defined in the center-of-mass frame to be 

£i = ^ ( l?_u l§) ■ (189) 
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In the case of the reactions (|180|) and (|181jl . i.e. with soft pion emission, we have to 
replace (|184|) with the three-particle phase space 

d 5 $(g + p; g',p', k) = ^l^l^( 2 ^S(q +p + k — q'-p'), (190) 

^% ^Po 

so that we obtain a differential cross section 



d 6^ir) = \M^^{ q +p-q\p\k) 

L 2(s-M 2 ) v /A(s,M 2 ,-Q 2 ) 
defined analogously to equation (|183j) . The amplitude is 

M? w) = e(N(p')it(k) n + (q')\e L ■ J|p(p)>, (192) 

and the superscript Nn labels the two final state possibilities mr or p7r~. Integrating out 
the angular dependence of the soft pion and the invariant mass W 2 of the final nucleon-pion 
system up to some (not too large) value W^., we further define 

w 2 

A2 (Nit) T* r j5 (Mr) 

d cr L _ f A ^2 r .„ dV 



w dn„ ——±-— . (193) 



dA 2 / / " dd)dA 2 dW 2 dtt 



Here, the variable fl n denotes the solid angle of the soft pion in the center-of-mass frame 
of the final nucleon-pion system. If the soft pion momentum k is sufficiently small, we can 
approximate the amplitude M.^ by its threshold value, and perform the phase-space 
integration in ()193|) to obtain 



d<P dA 2 32tt 2 (s - M 2 )A(s, M 2 , -Q 2 ) ' 
with the phase-space function 



(194) 



O 



th 



(195) 



$(^ max ) = ±J-£^(W«- - ^ h ) 3/2 

5.2 Transverse spin asymmetry 

Besides the longitudinal cross section itself, there exist predictions for the so-called trans- 
verse spin asymmetry in 7r + production (|182|) with a polarized target proton [23 [2^. There 
were arguments brought forward that this observable is particularly useful, because it is 
less sensitive to higher-twist corrections and next-to-leading order corrections in the strong 
coupling a s [29]. The definition of this asymmetry is 
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where S± is the component of the proton's spin vector that is transverse to q in the 
center-of-mass frame. A splitting of the squared amplitude into spin-dependent and spin- 
independent parts yields 

EiM n) i 2 4 n W*o)+4 n W^)i£idW: (197) 

S> 

so that one finds 

2s (n) 

A a = — l w (198) 

It is now straightforward to define an appropriate generalization of this observable for 
the equivalent process with soft pion emission: 

ANn = _L / d0---- - / d0^— / d0---- = ^Lp,, (199) 
\S ± \\J "d0dA 2 J r d0dA 2 \7 ^d0dA 2 7rs (^)' 



where the functions Sq^ and arise from the following decomposition of the threshold 
amplitude: 

E iM^r « ^[4^(^, A 2 ) + s^W A 2 )|5J sin# (200) 

The constant of proportionality not written explicitely equals that one in (|197j) , it contains 
the distribution amplitude of the n + and the strong coupling a s , for example. Supple- 
menting the factor f 2 in the denominator serves to keep the functions s$ and 
dimensionless. 

Of course, from an experimental point of view, these asymmetries are useful only when 
the detection of the soft pion is guaranteed. If this is not the case, one should formulate 
an asymmetry with the cross sections of all three processes, (|180j) to l|182j) . added up. In 
terms of the functions Sq and S\, such an asymmetry reads 

- 4 »+"'°+-- 7 s (») +(s (».») +s r- , )*(^)//r 1 ' 

5.3 Amplitude at threshold 

Now we come to the amplitude M.^^ of hard pion production with soft pion emission. 
For the sake of generality, we give the amplitude for arbitrary pion isospins a and c: 

e(N(p')7c a (k)rc c ( q ')\e L - J\N(p)) 

+(c + 5 ac /3 + c-Ke^iH^ + H™Mn) 

+ {c + Be acd r d /3 + c- (5 ac r 3 - 5 a3 r c ) + H ( ~^Mn)}^ 5 U, (202) 
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where we abbreviate 

± - 1 1 

- x _ e + fio ± a; + e-fio- (203) 

Evaluating this expression for the isospin combinations needed in ir + production, and 
writing the amplitudes with the help of two functions A^^ and Cw n in the form 

i 

M K L = — — J du Y^ U W ) i A Nnn + J 75«(p), (204) 



we obtain 



i 

A '™° = -^hi^rm-T^m^™^ (205 » 

i 



for the n7r° final state, and 

i 



V- = fj^(7^-7^)^ x ^ + 6 ^ x ^ (207 » 



-1 
1 



<V- - %fte(j^-j^)\i®M<) + m*,tt)l (208) 

-1 

for the p7r~ final state. 

Squaring the amplitude and summing over the final nucleon spin S', we obtain the 
functions s$ and Sj that we have introduced within the decomposition of the squared 
amplitude (127)01) : 



4^ = 4[(Re^) 2 + (Im^) 2lM(1 ( " 



M + m n 



+ ^ + ™^ Re(A N7r C* Nn ) - jL [(ReC^ f + (ImCW) 2 ] (209) 

*r = -+£^MAM. (210) 

Note that when setting the pion mass equal to zero everywhere in the last expression, we 
recover the known formula for the reaction without soft pion emission. 
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5.4 Numerical results 



5.4.1 Ratio of the longitudinal cross sections 

The longitudinal cross section for usual hard ir + production to leading twist and leading 
order in QCD was given in references [3H[21|. Here, we present the ratios of the unpolarized 
cross sections 



d<rf *> _ T dVf 



dA 2 J "d0dA 



2 



(211) 



o 



and their corresponding counterpart without soft pion emission, respectively. These ratios 
must be calculated for a fixed direction of q', which implies the following relation between 
two- and three-particle final state variables A 2 and to: 



x B (Wl -M 2 + t )-t (± 
1-xb W 5 



A 2 = - * v " th - J V ^ + [^-). (212) 



With this specification of A 2 , we can define the ratios as functions of xb and to'- 

R»,fei.,HU-^r^j " 4») (lB , to) n— (213) 

Concerning the models for the GPDs that enter the calculations of and s n \ we 

refer to the review [lj; relevant formulas are summarized in appendix El For the pion 
distribution amplitude, throughout we use the asymptotic form 

<PM = <p as {u) = - 4 {l-u 2 ). (214) 

Figure El shows the ratio as a function of xb for three different values of t . We 
find that for a large range of t , the soft-pion contamination is roughly 10% for each 
individual channel. In situations where we have to add up both contributions, we arrive 
at a soft pion contamination that accounts for 15% to 30% within the presented range of 
to- The J¥ max -dependence near threshold can easily be deduced from our results given at 
W max = 1.15 GeV 2 , because it is incorporated explicitely in (|213|) through the phase space 
function $(W max 



5.4.2 Transverse spin asymmetry 

The results for the transverse spin asymmetries of the processes with soft pion emission 
(jl80j) and (|181jl are shown in figure The values for these asymmetries turn out to be 
very small compared to those of the pure process (|182L see figure |H1 The small size of the 
asymmetries ^4^^ can be traced back to a particularly small value of sf^^ that arises from 
a significant cancellation within the following terms: 

oc ImAjsfnReCNTT — ReAAr^ImCArTr- (215) 
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Figure 6: The two ratios of the unpolarized longitudinal cross sections with and without 
soft pion emission, R n7T o and -R p7r -, as functions of xb for the values t = —1 GeV 2 , t = 
-0.3 GeV 2 and t = -0.02 GeV 2 . Lower curves belong to smaller values of —t . The 
photon virtuality is Q 2 = 10 GeV 2 and W max = 1.15 GeV. The curves are plotted up to the 
maximally allowed value of the Bjorken variable, which is XB max = 2/(1 + a/1 — 4M 2 /t ). 
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Figure 7: The transverse spin asymmetries of the processes p + 7* — > n + 7Tg oft + 7r + and 
p + 7* — > p + 7r~ ft + 7r + for the momentum transfers to = —0.3 GeV 2 (lower curves) and 
to = — 1 GeV 2 (upper curves). The photon virtuality is chosen as Q 2 = 10 GeV 2 . 



In contrast, such a cancellation does not happen in the calculation of s\ , because in the 
pion-pole model for E, (see appendix EJ the imaginary part of C n vanishes, whereas A n 
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Figure 8: The transverse spin asymmetry for the pure hard n + production, A n (dashed 
curve), and for the hard ir + production with soft pion admixture, ^4 n +n7r + P 7r- (solid curve). 
The values of the momentum transfer are to = -0.02 GeV 2 (left), t = -0.3 GeV 2 (middle) 
and to — — 1 GeV 2 (right). In the case of A n+n7T o +p7T - , the maximal invariant mass is chosen 
to be W max = 1.15 GeV. The photon virtuality is Q 2 = 10 GeV 2 . 



has non-negligible real and imaginary parts, 

-1 
-1 

Moreover, the values of s N ^ are larger than those of roughly speaking by a factor 
two or three. Since they occur in the denominator, An-k = ^s^^ /(irs^^), this gives a 
further reduction of the asymmetry An-k- 

Let us now come to the asymmetry A n+n7T o +pn - of the combined processes, see equation 
()20ip . Our discussion of the functions s and si above already indicates that the inclusion 
of soft pions leads to an certain reduction of the asymmetry. This result is shown in figure 
IH1 For an invariant mass integrated up to W max = 1.15 GeV, we find a downwards shift of 
the curves of about 10%, respectively. 

6 Summary 

In the following, we summarize our statements about soft pion emission in hard exclusive 
reactions: 
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• We have formulated a parametrization of twist-2 lightcone operators between initial 
nucleon and final nucleon-pion state. It has turned out that four functions, which we 
call irN partem distributions, are needed for each of these matrix elements. Further, 
we have shown how this number is reduced to two when the pion is exactly at the 
production threshold. 

• Next, we demonstrated on several examples that the moments of the ttN distribu- 
tions are polynomials in the skewedness variables £ and a. The coefficients in these 
polynomials are essentially the form factors (or, following the terminology of elec- 
troweak pion production, the "invariant amplitudes") of pion emission induced by 
corresponding local twist-2 operators. 

• We have derived soft-pion theorems for the irN distributions which represent the 
leading terms in an expansion in small pion momentum and mass at given momentum 
transfer. The nontrivial ingredients of the final expressions are nucleon and pion 
GPDs as well as the pion distribution amplitude. If the momentum transfer is fixed 
to be large compared to the pion mass, we find that our results agree with Guichon 
et. al. |Hj. In the opposite case, i.e. when the momentum transfer is small, we have 
argued on the example of certain moments that our results are consistent with the 
leading order of ChPT. 

• Moreover, we have provided explicitely a parametrization and the soft-pion theorems 
for the resulting form factors of pion emission from the nucleon induced by the energy- 
momentum tensor. 

• Finally, we have given analytical results for the amplitude of hard pion production 
with soft pion emission and some numerical estimates for the particular case of hard 
7r + production using certain GPD models. In the presented kinematical region (xb > 
0.05, 0.02 < -Ag/GeV 2 < 1, VF max = 1.15 GeV, Q 2 = 10 GeV 2 ), we have obtained 
that the contamination of the longitudinal cross section of the reaction 7* + p — > 
7r + + n caused by soft pions amounts to roughly 10% within each of the two possible 
channels 7* + p — > n + + n + 7Tg oft and 7* + p — > n + + p + 7T~ ft . The effect of the 
soft-pion channels on the transverse spin asymmetry is a downwards shift of roughly 
10%. Further, we have observed that the transverse spin asymmetry of the individual 
soft-pion reactions is about one order of magnitude smaller than that one which is 
predicted for the familiar process without soft pion. 

We conclude that soft-pion theorems can provide useful estimates of the soft-pion contam- 
ination in hard exclusive reactions. The numerical studies for the contamination of hard 
7r + production indicate that for an accurate interpretation of experiments an appropriate 
separation of the soft-pion channels is required. Second, we point out that in principle, 
according to the soft-pion theorems, hard exclusive reactions with soft pion emission can 
serve as an additional source to extract information about GPDs. Therefore, it is worth 
to consider corresponding experiments as well as further theoretical studies. 
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Appendix 



A GPD models 

In this appendix, we summarize the models that we have used for the isovector nucleon 
GPDs which enter the presented calculations of hard ir + production with soft pion emis- 
sion. They are taken from reference pQ, to which we refer for further explanations and 
motivations. 



A.l Isovector GPD 

The isovector and isoscalar GPDs and are linear combinations of the up and 

down quark GPDs in the proton, 



H (V) =R u_ H d H (S) =R u + R d_ 



(218) 



The models for the quark GPDs H u and H d contain the so-called double distribution term 
and the .D-term. The latter is supposed to be flavor-independent, hence it cancels in the 
isovector combination. Further, one assumes a factorized ansatz for the quark GPDs H q , 
q = u,d, so that 



HM{x,t,t) = Hl D {x,t) 



F? /P (t) 



H d DD {x^)Ff /p {t) 



where the double distribution representation is 



(219) 



dp J daS(x- (3- of) h((3, a) H q ((3, 0,0). 

-1 -(1-|/3|) 

Here, the profile function h is chosen as 



h(/3, a) 



3(1- 
4 (1-|/3|)3 



a 



(220) 



(221) 



and H q (x, 0,0) denotes the forward limit that is related to the usual quark distribution 
functions q(x) and antiquark distribution functions q(x) according to 



H q (x, 0, 0) = q(x) 9(x) - q(-x) 6(-x). 



(222) 



For the numerical calculations, we use the leading order parametrization MRST2001LO 
from reference |32]- The second ingredient of the factorized ansatz for H is the Dirac form 
factor . It is defined according to 



(N f (p')\^ q ^ q \N f (p)) = U' 



2M 



u. 



f =P,n, 



(223) 
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and related to the nucleon Sachs form factors via 



W<« - |[J?"(f) + F? /f (t)\ - \[F? f (t) + F 2 d/f (t)} 



G { i\t) 



F? /f (t) + 



tF^ /f {t) 
AM 2 



F? /f (t) + 



tFj /f {t) 
AM 2 



(224) 
(225) 



As input for numerical calculations, we have used the empirical fits of Brash et. al. |33j for 
the proton Sachs form factors and those of Bosted [34J for the neutron ones. 



A. 2 Isovector GPD 

For E^ v \ we take the very simple model 

E( y \x,^t) = [El D {x,i)-E d DD {x,i)]G D {t) 
with the double distribution representation 

1 1-1/31 

E q DD (x, = J d/3 J da5(x-(3- aO E*(/3, 0, 0) h{/3, a). 



(226) 



(227) 



-1 -(1-1/31) 

Concerning the forward limit of E q , the ansatz 



K 



with 



E u (x, 0, 0) = —u v (x) 9(x), E d (x, 0, 0) = K d d v (x) 9{x) 



k u = 2k p + K n K d = k p + 2n n 



(228) 



is assumed. Here, k v = 1.793 and n n = —1.913 are the anomalous magentic moments of 
the proton and neutron, and q v (x) denotes the valence quark distribution q(x) — q(x). In 
the t-dependent factor of the ansatz for E^ v \ Gd is the dipole form factor 



G D (t) 



(l-i/0.71 GeV 2 ) 2 ' 



(229) 



A. 3 Isovector GPD 

The GPD = H u — H d is modeled by the double distribution ansatz 

H«(x,{i,t) = H q DD (x,{i)Fl(t) 

with 

1 1-1/31 

H q DD {x^) = y d/3 J da6(x-p-a£)h(p,a) Aq v (J3), 

-1 -(1-1/31) 



(230) 



(231) 
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where Aq v = [Aq(x) — Aq(x)]9(x) is the polarized valence quark distribution. For the nu- 
merical input of Aq v we have used the leading order analysis LSS2001LO given in reference 
[35] . Finally, the axial form factors are approximated by the dipole form of F^\t), 

with an axial mass Ma ~ 1 GeV. 

A.4 Isovector GPD 

The GPD EM is modeled by the pion pole form 

= G^\t)-^(x/0 e(t - \x\), (233) 
where the pseudoscalar form factor is assumed to be 



ml — t 



G { p'(t) = " ' V ' (234) 



and for the pion distribution amplitude we use its asymptotic form as already given 
in K0 



B Soft-pion theorems for non-threshold pion-nucleon dis- 
tributions 

In section H31 we have shown the soft-pion theorems for the threshold tiN distributions. 
Here, we give the corresponding general results, i.e. when the pion slightly deviates from 
the threshold. Throughout the following expressions, we shall imply the relations 

u = 2M 2 + ml - W 2 - t + A 2 , (235) 



a = a(l-f), (236) 
(P + p ) ■ n 2 — a 2 

and 

(p — p' + k) ■ n , 
p t -n= 2 = £ + «, (238) 

which make the dependence of the ttN distributions Hf'^ on the set of variables x, £, A 2 , 
a, t, and W 2 explicit. 
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B.l Isovector vector ttN distributions 



H[ +) = E^\x^^) + ^%[H^\x^^)+E^\x^^)] (239) 

u — M z 

H^ +) = (240) 



H[-> = |°[BW( I , & , t )-0 I .&. t )l-^ ! ^[« (v) (^fA J ) + ^( I ,|,A 2 )] 

+ ^~M V) (J?-, J-, o) 6( Pt .n-\x\) (243) 
"4 - * \Pt-n p t -n J 

Ht ] = - H(V) ( X > *>> + gOOfo ^ A 2 ) + E<V(x, A 2 ) 



M 2 M 2 



(244) 



- -(ra^-^J^-^ 2 ' (245 » 

B.2 Isoscalar vector and gluon ttN distributions 

H (o,G) = ^, G ) ( ^ )A 2 ) + ^ . [H (s, G) ( ^ A 2 ) + ^ G )(x,e,A 2 )] 

u — M z 

+-^W >G) M", z^-, o) £fe • n - x) (247) 
m*-t \Pt-n pt-n J 

H^ G) = (248) 

^ = -(F^ + ^) |F<0 ' G,fe ^ A2) + E< °' G,(l ' f - A2)1 <250) 
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B.3 Axial vector tiN distributions 





= ^ v \x^,A 2 ) 


2M 2 a 
u-M 2 


#(0,+) 


= 




fr(0,+) 


/ M 2 


M 2 \ 


\W 2 -M 2 


u-M 2 ) 




} M 2 


M 2 \ 


\ W 2 — M 2 


u-M 2 ) 



H {sy \x,^A 2 ) (251) 

(252) 

Z& s > v \x,Z,A 2 ) (253) 
H (s ' v \x,Z,A 2 ), (254) 



^ = W^+Mj.BOOM,*) (255) 

^ - f 1 i M2 I M2 ^ CE (v) (x C A 2 ) 
Hs ~ \2 + W 2 -M 2 + u-M 2 )^ ^' A ^ 

2M 2 fr\ 
+ , , ^ | W - H) (257) 



M 2 M 2 



& - Hra + ^r ( ^ A2) - (258) 

C Soft-pion theorems for the form factors of pion emis- 
sion induced by the local vector current 



In equation l|3*%|) of section 13.31 in the main text, we have parametrized the matrix element 
for pion emission induced by the local vector current in terms of form factors A^. As 
explained in this section, the soft-pion theorems for these form factors can be obtained 
from the first moments of the rcN distributions functions that have been given in 

appendix [Bj In this way, we obtain for the isoscalar and isovector even form factors: 

8 

Y^U' A?' +) T?U 



u 



M 2 



1 + " 



II -' - M 2 u-M 2 
' M 2 M 2 



W 2 -M 2 u-M 2 



[F} sy) (A 2 ) + F^ v \A 2 )]kY \ IsU + 0(e). (259) 
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To the given accuracy, current conservation is fulfilled: 

i 

For the isovector odd form factors, we obtain 

u' 



(260) 



i=i 



2g, 



+ [-FZ\t) + F^>(^) + F^>(^)}M r 



A^ + 



2M 2 [F 1 (y) (A 2 ) + Ff ) (A 2 ^ 



M 2 



+ 



AM 2 
m 2 —t 



M 



M 



W 2 -M 2 u-M 2 



M 



M 



II' 2 - M 2 u-M 2 
M 2 M 2 



W 2 -M 2 



— M 2 



F 2 (y) (A 2 )^ 

[FP(A 2 ) + FM(A 2 )]hfy i r,U + OL~) 



(261) 



Current conservation at small momentum transfer is immediately fulfilled, but for the mod- 
erate momentum transfer, the pseudoscalar form factor Gp cannot simply be approximated 
by its pion pole form, and so we arrive at 



M 



G A (t) 2M + G P {t) 



t 



2M 



U' l5 U + 0(e) (-t>£ 2 ). (262) 



But since the remaining combination of axial and pseudoscalar form factor is proportional 

to m 2 , 



G A {t) 2M + G P {t) 



t 



2M 



7sy£/ = (N(p')\d-A a \N(p)) = f n ml(N(p')\<S> a \N(p)) (263) 



and pion-pole enhancement at such large t is excluded, it is reasonable to assume GU(£) 2M+ 
Gp(t) tj (2M) = 0{e). Hence for any momentum transfer in the considered region —t < M 2 
we arrive at 

A" U'A\' ] TfU = 0(e), (264) 



so that finally, current conservation is fulfilled. 
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